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We present the first NLO QCD analysis of single-diffractive dijet production in proton-antiproton 
collisions. By comparing the ratio of single- and non-diffractive cross sections to data from 
the Tevatron, the rapidity-gap survival probability is determined as a function of the momentum 
fraction of the parton in the antiproton. Assuming Regge factorization, this probability can be 
interpreted as a suppression factor for the diffractive structure function measured in deep-inelastic 
scattering at HERA. In contrast to the observations for photoproduction, the suppression factor 
in proton-antiproton colhsions depends on the momentum fraction of the parton in the Pomeron 
even at NLO. 
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1. Introduction 

In high-energy hadronic coUisions, a large fraction of the events is produced diffractively and 
contains one or more leading (anti-)protons or excited hadronic states Y with relatively low mass 
My, which are separated from the central hard process by rapidity gaps. Hard diffraction has been 
studied intensely in deep-inelastic ep scattering (DIS) at HERA, where it can be understood due 
to the presence of the large squared photon virtuality Q^ in terms of the Pomeron structure func- 
tion F^{x,Q^,^,t) = Y^iffip ® ^y*i' factorizing into diffractive parton density functions (DPDFs) 
ffj and perturbatively calculable partonic cross sections o^i [jl]]- Assuming in addition Regge 
factorization, the DPDFs f^, = fi/p{p = x/^,Q^)fpip{^,t) can be parameterized in terms of a 
Pomeron flux factor fp/p, which depends on the longitudinal momentum fraction of the pomeron 
in the proton ^ and the squared four-momentum transfer at the proton vertex t, and PDFs in the 
Pomeron fj/p, which depend on the momentum fraction of the partons in the pomeron j8 and an 
unphysical factorization scale Q^, and fitted to diffractive DIS data [§]. Diffractive dijet production 
data add further constraints, in particular on the gluon distribution [^]. At high <§, a subleading 
Reggeon contribution must be taken into account. 

As Q^ decreases, the applicability of perturbative QCD must be ensured by sufficiently large 
jet transverse energies Eti.2- From the absorption of collinear initial-state singularities at next-to- 
leading order (NLO) [Q], the photon develops a hadronic structure and contributes not only with 
direct, but also resolved processes ^. These processes occur also in hadron-hadron scattering, 
where factorization is broken due to the presence of soft interactions in both the initial and final 
states [^, and factorization breaking is indeed observed in dijet photoproduction at NLO [^. Sim- 
ilar effects may also arise in dijet production with a leading neutron ^. Soft rescattering must 
be understood, if one wishes to exploit processes like diffractive Higgs production, which offers a 
promising alternative to the difficult inclusive diphoton signal in the low-mass region [P]. In this 



article, we present the first NLO analysis of diffractive dijet production in hadron collisions QlOl]. 



2. Diffractive dijet production at the Tevatron 

Dijet production with a leading antiproton has been measured by the CDF collaboration at 
the Tevatron in the range |f| < 1 GeV^, used also by the HI collaboration for the extraction of 
their PDFs 2006 Fit A,B and 2007 Fit Jets, but at shghtly larger ^ e [0.035; 0.095] ^. From the 
ratio R{xp = j8(^) of single-diffractive (SD) to non-diffractive (ND) cross sections as a function of 
Xp = L(^rr/v^^ '^'» integrated over En > 1 GeV and |t],| < 4.2, an effective diffractive structure 
function Ffj{l5)=Rx F^/^ was extracted assuming /^/^(.x, Q^) = x[fg/p + 4/91,1 fq_/p] with GRV 
98 LO proton PDFs and compared to older HI diffractive structure functions at fixed Q^ = 75 
GeV^. These were found to overestimate the measured diffractive structure function by about an 
order of magnitude, while the correction to be applied to the HI measurements with My < 1.6 GeV, 
when compared to the CDF measurements with leading antiprotons, should only be about 23% [Q]. 

In a subsequent publication, the CDF collaboration analyzed in a similar way data at ^/s = 
630 GeV, not only 1800 GeV, in order to test factorization, Pomeron flux renormalization, and 



rapidity-gap survival models [12], but found no significant energy dependence [13]. In addition 
they restricted the |f|-range to below 0.2 GeV^ and the average transverse jet energy to Et > 10 



NLO QCD analysis of single-diffractive dijet production at the Tevatron 



Michael Klasen 



10 



>0^ -2 

10 



10 



pp -^ (p+)2jets+X at Vs = 1800 GeV 

0.035 < ^ < 0.095 



h 






***** 

h-«-H CDF PRL 84 
w^M "NLO 
LO 

.NLO HI 2006 Fit A 

^^ ,NLQHl,20Q7Fit.fets ,^^^ ,^^^ ,^^^ , /<- 
3.5-3.25 -3 -2.75-2.5-2.25 -2 -1.75-1.5-1.25 -1 



0.3 ^ 
0.25 - 

g-0.2 - 

z 
■To, 

X 

a. 

?oi 0.1 - 



pp -> (p-l-)2jets-HX at Vs = 1800 GeV 



T 



0.035 < ^ < 0.095' 



if 



**r 



0.05 







iH-^ CDF PRL 84 / NLO [ft^ * 

|-^>H CDF PRL 84 / LO >»»<* 

:i-n CDF PRL 84 / NLO HI 2006 Fit A 
- , Cp^PR|:. 84/NLpHl, 2007 Fit ^ets , 
3.5-3.25 -3 -2.75-2.5-2.25 -2 -1.75-1.5-1.25 
logio(\) 



Figure 1: Left: Ratio R of SD to ND dijet cross sections as a function of the momentum fraction of the parton 
in the antiproton, computed at NLO (with three different DPDFs) and at LO and compared to the Tevatron 
Run I data from the CDF collaboration. Right: Double ratio of experimental over theoretical values of R, 
equivalent to the factorization-breaking suppression factor required for an accurate theoretical description of 
the data (color online). 

GeV, which has the advantage of removing the infrared sensitivity of two equal cuts on £71,2 [0]- 
In both publications, the extraction of Fjj was performed on the basis of leading order (LO) QCD 
and assuming that the f-channel gluon exchange in the partonic cross sections dominates. 



3. NLO QCD analysis 

Our NLO QCD analysis is based on previous work on inclusive dijet photoproduction [[T^], 
where the resolved photon contribution can be directly applied to dijet hadroproduction [|l6|]. It 
takes into account the convolution of Pomeron flux factors and PDFs, modem CTEQ6.6M proton 
and HI 2006 and 2007 Pomeron PDFs with Q^ = Ej^ on an event-by-event basis, and the complete 
set of parton-parton scattering cross sections. Jets are defined by a cone with radius R = 0.7 as in 
the CDF experiment and imposing a parton separation of /?sep = 1-37?. The normalized control 
distributions \/odo /dEj and l/ada/df] (not shown) are in good agreement with the SD and ND 
measurements, in particular for Ej2 > 6.6 GeV removing the infrared sensitivity in the first CDF 
analysis and for y^ = 630 GeV in the second CDF analysis [ 10 1 . However, the ratio R of SD to ND 
dijet cross sections as a function of Xp shown in Fig. [I] (left) for y/s = 1800 TeV and |?| < 1 GeV^ 
significantly overestimates the data for all employed diffractive PDFs and even more at NLO than at 
LO, due to an average ratio of SD to ND /T-factors of 1.35 (1.6 at ^/s = 630 GeV). Remember that 
all theoretical predictions should be divided by a factor of 1 .23 to correct for diffractive dissociation 
included in the HI DPDFs. 

The double ratios R^^^ /R^^)^'^ in Fig. |T] (right) and Fig. ^ can be interpreted as suppression 
factors or rapidity-gap survival probabilities, which are smaller than one due to soft rescattering. 
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Figure 2: Double ratios of experimental over theoretical values of R, equivalent to the factorization-breaking 
suppression factor required for an accurate theoretical description of the data from the Tevatron at ^/s = 630 
(left) and 1800 GeV (right) (color online). 



Their qualitative behavior in these three figures is very similar: we observe a dependence on Xp with 



a minimum at \ogyQ{xp 



4.5 (.Xp ~ 



0.032), a rise towards smaller xp by up to a factor of five, 



a smaller rise towards larger xp, and an appreciable dependence of the suppression factor on the 
chosen diffractive PDFs. The corresponding suppression factors for F^(j8) (not shown), computed 
with the same assumptions on Ffj^ as those taken by CDF, are minimal with the value ~ 0.05 at 
j8 = 0.5 and rise to ~ 0. 12 for ^^ = 1800 GeV and to ~ 0. 1 for ^/s = 630 GeV at j8 = 0. 1 and with 
HI 2006 Fit B. Here, the effect of NLO corrections is partially compensated by the simplifications 
inherent in the calculation of Fjj. The strong rise of the suppression factors to values of ~ 0.25 
for j8 < 0.05 [logiQ^xp) < -2.5] at ^^ = 1800 GeV may not be conclusive, since the HI DPDFs 
are largely unconstrained in this region. On the other hand, the rise at log]o(.Xp) > — 1.5 can be 
attributed to the Reggeon contribution, that should be added at large ^ . Note also that the variation 
of the suppression factor is considerably reduced for the HI 2007 Fit Jets, which has been obtained 
with additional constraints from diffractive DIS dijet production. 

The extraction of F/j(j3) from R{xp) is based on the assumption that the latter is only weakly 
(^-dependent and can be evaluated at an average value oi ^ = 0.0631. This weak (^-dependence 
is indeed observed in the newer CDF data and also in our theoretical calculations, which reflect 
the (^-dependence of the HI fits to the Pomeron flux factors fp/p{^,t) <=< (^^'" with m ~ 1.1 (0.9 
in the CDF fit to their data). The (small) difference of the theoretical (1.1) and experimental (0.9) 
values of m can be explained by a subleading Reggeon contribution, which has not been included 
in our predictions. To study its importance, we have computed the ratio of the Reggeon over the 
Pomeron contribution to the LO single-diffractive cross section at ^/s = 1800 GeV. The Reggeon 
flux factor was obtained from HI 2006 Fit B and convolved, as it was done in this fit, with parton 



densities in the pion [17]. Very similar results were obtained for the HI 2007 Fit Jets Reggeon flux. 
On average, the Reggeon adds a 5% contribution to the single-diffractive cross section, which is 



NLO QCD analysis of single-diffractive dijet production at the Tevatron Michael Klasen 

smaller at small xp = ^{3 and <§ (2.5%) than at large xp and £, (8%). This corresponds to the graphs 
shown in Figs. 5 (i^ = 0.01) and 6(^= 0.03) of Ref. [|], e.g. at Q^ = 90 GeV^. While the Reggeon 
contribution thus increases the diffractive cross section and reduces the suppression factor at large 
Xp, making the latter more constant, the same is less true at small values of ;cp. 

Any model calculation of the suppression factor or rapidity-gap survival probability must try 
to explain two points, first the amount of suppression, which is ~ 0.1 at j8 = 0.1, and second its 
dependence on the variable j8 (or xp). Such a calculation has been performed by Kaidalov et al. 
[ir2|]. In this calculation, the hard scattering cross section for the diffractive production of dijets 
was supplemented by screening or absorptive corrections on the basis of eikonal corrections in im- 
pact parameter space. The parameters of the eikonal were obtained from a two-channel description 
of high-energy inelastic diffraction. The exponentiation of the eikonal stands for the exchange of 
multi-Pomeron contributions, which violate Regge and QCD factorization and modify the predic- 
tions based on single Pomeron and/or Regge exchange. The obtained suppression factor is not 
universal, but depends on the details of the hard subprocess as well as on the kinematic configu- 
rations. The first important observation from this calculation is that in the Tevatron dijet analysis 
the mass squared of the produced dijet system Mjj = Xpp^s as well as t, are almost constant, so 
that small j8 implies large Xp. The second important ingredient in this calculation is the assumption 
that the absorption cross section of the valence and the sea components, where the latter includes 
the gluon, of the incoming proton are different, in particular, that the valence and sea components 
correspond to smaller and larger absorption. For large Xp or small j8 , the valence quark contribu- 
tion dominates, which produces smaller absorptive cross sections as compared to the sea quark and 
gluon contributions, which dominate at small Xp. Hence the survival probability increases as Xp in- 
creases and j3 decreases. The convolution of the j8 -dependent absorption corrections with older HI 



DPDFs 1 18] led to a prediction for F^{P), which was in very good agreement with the correspond- 
ing experimental distribution (see Fig. 4 in [|l2|]). A similar correction for soft rescattering of our 
single-diffractive NLO cross sections based on the more recent DPDFs of HI should lead to a very 
similar result. An alternative model for the calculation of the suppression factor was developed by 



Gotsman et al. []19|]. However, these authors did not convolve their suppression mechanism with 
the hard scattering cross section. Therefore a direct comparison to the CDF data is not possible. 

At variance with the above discussion of diffractive dijet production in hadron-hadron scat- 
tering, the survival probability in diffractive dijet photoproduction was found to be larger (~ 0.5 
for global suppression, ~ 0.3 for resolved photon suppression only) and fairly independent of j8 
1^, ^^. This can be explained by the fact that the HERA analyses are restricted to large values of 
X'y>0.l (as opposed to small and intermediate values of Xp = 0.02 ... 0.2 at the Tevatron), where 
direct photons or their fluctuations into perturbative or vector meson-like valence quarks dominate. 
The larger suppression factor in photoproduction corresponds also to the smaller center-of-mass 
energy available at HERA. 
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